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1 . INTRODUCTION 


One  of  the  problems  encountered  in  the  study  of  impact  phenomena  is 
the  determination  of  a velocity  below  which  an  object  will  fail  to 
perforate  a barrier  or  some  type  of  protective  device.  This  determina- 
tion is  of  prime  importance  in  the  design  of  protective  structures  in 
nuclear  power  plants,  for  the  development  of  safety  features  in  the 
automotive  industry,  for  evaluating  the  effectiveness  of  military 
vehicle  armor  and  in  any  problem  area  where  an  impact  can  cause  damage. 
This  velocity  is  commonly  referred  to  as  the  critical  impact  velocity. 

The  techniques  available  to  determine  this  velocity  can  be  classed 
as  either  deterministic  or  probabilisti*..  In  the  former  category,  a 
-Umit  velocity  (denoted  as  Vj^)  is  determined  from  physical  principles 

(the  conservation  laws  and  material  constitutive  relations)  but  because 
of  the  complexity  of  the  governing  partial  differential  equations, 
simplifications  arc  introduced  which  generally  require  empirical 
determination  of  one  or  two  constants.  In  the  probabilistic  approach, 
models  are  built  relying  on  a substantial  base  of  data  consisting  of 
the  object's  striking  velocity  and  either  its  residual  velocity  or  a 
statement  of  either  defeat  or  non-defeat  of  the  barrier.  The  resulting 
critical  velocity  is  most  commonly  expressed  as  a i.e.,  a striking 

velocity  for  which  there  exists  a 50%  probability  of  perforation  of 
the  barrier.  In  determinations,  a statistical  approach  is  employed 

wherein  the  response  is  quantal  and  a sensitivity  test  can  be  applied 
to  the  data.  The  simplest  approach  requires  determination  of  a mean 
velocity  for  three  data  points  above  and  below  the  estimated  critical 
impact  velocity. 

In  practice,  time  and  economic  constraints  Umit  the  quantity  of 
data  obtained  so  that  there  is  always  less  than  that  desired  by 
either  deterministic  or  probabilistic  models  for  reliable  estimation 
of  critical  impact  velocities.  After  a brief  discussion  of  various 
models  currently  in  vogue  and  their  advantages  and  limitations,  this 
paper  discusses  procedures  designed  to  extract  reliable  estimates  of 

or  from  sparse  data  sets. 

To  associate  the  limit  velocity  with  the  technique  used  to  obtain 
the  critical  impact  velocity  this  report  defines  acronyms  as  follows: 

VL  a limit  velocity  obtained  from  a graphical  determination 
or  an  empirical  formulation. 

F’BL  a protection  ballistic  limit  determined  from  a special 
experimental  technique. 

VLL  a limit  velocity  based  on  an  analytical  representation 
of  striking  and  residual  velocity  data  developed  by 
Mr.  J.  P.  Lambert  of  the  BRL. 
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' tMllist.v  limit  Vilocify  ui*c-d  when  it  is  miesMiry  In 

.upiMiMt  tor  mixed  rcnponixc  hi  the  dat:i. 


11.  I.IMII  VI.UK:m',  V,  I urn:  KM  I NAT  I ON 

M 

ihe  liuit  velocity,  V is  detorminod  by  observint  The  rclutiunship 
iicttaccn  the  striking;  velocity,  of  the  projectile  und  its  residual 

velocity,  V^,  for  several  striking  velocity  levels  above  an  estisiatcd 

critical  velocity.  Starting  with  a very  high  velocity  the  test  is 
repeated  until  a very  low  residual  velocity  Is  observed.  A plot  is 
:j.ide  of  the  residual  velocities,  V^,  as  a function  of  the  striking 

velocities,  \\..  by  drawing  an  “cyobair*  curve  through  the  data  points 

and  extrajiolat  ing  this  curve  to  Vj^  ■ 0,  a value  for  V'g  Is  found  below 

rtlilcli  the  projectile  will  fail  to  defeat  the  target.  This  striking 
velocity  IS  defined  as  the  limit  velocity,  Vj^,  Figure  1 is  a graphical 

represent  at  ion  of  such  a curve.  This  methoi!  is  straightforward  and 
involve*,  no  statistical  procedure.  With  u careful  selection  of  striking 
velocities  the  can  be  determined  from  very  few  tests.  However,  if 

"iixed  rc-ults  occur  as  the  striking  velocity  is  reduced,  then  another 
'icihiHl  must  he  used  to  obtain  the  limit  velocity. 


m.  rpoTi.CTiON  BAi.i.isric  i.imit,  i’bi. 

1)10  jirotcclion  ballistic  limit,  PHI.*,  adds  a protection  critviioit 
to  the  evaluation  of  the  crlticol  impact  velocity.  A comjiletc 
poneLraiion  is  obtained  when  a fragment  or  fragments  of  either  the 
impacting  projectile  or  armor  arc  caused  to  be  thrown  from  the  plate 
with  .sufficient  energy  to  perforate  a .sheet  of  .50mm  (.020  in.)  thick 
iiural  placed  parallel  to,  und  15cm  (6  In.)  to  30cm  (12  in.)  behind 
the  test  Item.  All  other  penetrations  arc  partial  penetrations.  Since 
in  actual  experimental  tests  variation  in  results  do  occur,  further 
cvaiuatiun  of  the  type  previously  described  mu.st  be  employed. 


* . ...  iru’ii'iK,  'TiJ>u>r.  rat.  Con  oj' /ir^r  Oj  tHeoL  and  Hi^h  Danait^- 

r.iLlCtttlo  siraoxiron  .^d>oratorico  >^’norafu:t4n  Hepor:  .V. 
U-i,  l.ri,  (AO  ■518.39.11.) 

• o'ui  i\rT>utniiion  fuatiut/”*  irroof  Manual,  Vol.  I, 

/it  r'r : '.O'.  Prooing  Uixjund,  MD,  Januavj  IdtiO. 
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Figure  I,  Graphical  l^eterai nation  of  Linit  Velocity  fnm  Vj^/Vg  Data. 


IV.  LAMBLRT  LIMIT  VBLOCI'fY, 

The  analytical  representation  of  the  dependence  of  residual  velocity 
on  striking  velocity  for  a specified  projectile-target  configuration 
by  the  Lambert  Method^  is  of  the  form: 


A (V  P - V,  P'^^^P 


'LL 


* 0 1 Vg  < V, 
, Vg  > V 


where  the  imposed  constraints  are: 

p > 1 and  0 _<  A ^ 1 

The  parameters  A,  p,  and  are  determined  by  a non-linear  least  square 

regression  analysis,  A is  the  asymptote  to  the  Vj^-Vg  data  specified  by 

the  equation  = A(Vg).  The  parameter  p is  the  steepness  factor  whose 

value  influences  the  curvature  near  the  V » 0 section.  The  ballistic 

K 

limit  velocity  VLL  is  obtained  as  a by-product  of  the  fitting  technique. 
Also,  associated  with  the  derived  fit  of  the  form  to  the  data  is  the 
root  mean  square  error,  S,  or  the  standard  error.  The  advantage  of 
tliis  technique  is  that  it  uses  the  entire  data  whether  it  has  or 
has  not  a mixed  response. 


V.  BALLISTIC  LIMIT  VELOCITY  CALCULATIONS 

In  experimental  tests  where  the  only  observable  result  as  a function 
of  the  striking  velocity  is  the  perforation  of  the  target  or  the  lack 
thereof,  a statistical  approach  must  be  utilized  to  determine  the 
ballistic  limit  velocity.  In  statistical  literature  this  type  of 
quantal  data  (yes  - no,  success  - failure,  go  - no  go,  etc.)  is  denoted 
as  "sensitivity  data"  where  the  striking  velocity  is  the  stimulus  and 
the  result  is  the  response.  In  the  analysis  of  such  data  a normal 
integrated  response  function  is  used  as  a model  to  determine  the  ballis- 
tic limit  velocity;  that  is,  the  critical  impact  velocity  corresponding 
to  a fifty  percent  probability  of  success  in  defeating  the  target. 

In  some  statistical  methods  like  the  Logit**  and  Probit^  methods  the 
striking  velocities  and  the  number  of  rounds  at  each  striking  velocity 
are  predetermined.  In  other  methods  the  choices  are  made  sequentially 
as  the  experiment  progresses. 


V,  ?,  Lambert t G,  ri.  Jonaa,  "Towards  Standardization  in  Terminal 
dalliatios  Testing:  Velocity  Representation: ^ Ballisoia  Reaearoh 
uworatories  Report  No.  IdSl,  January  1976. 

^L).  J,  Finney,  "Tne  Estimation  of  the  EO  60  for  a Logiatia  Reaponae 
Curve",  Sanknya,  Vol.  12,  Parts  I and  2,  1962. 

^0.  J.  Finney,  "Probit  Analysis-A  Statistical  Treatment  of  the 
Sigmord  Response  Curve",  Cambridge  University  Press,  1962. 
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A.;  Moan  l.imit  Velocity,  Calculation 


A simple  and  straightforward  approach^  is  to  extract  from  a series 
of  experimental  results  the  three  highest  striking  velocities  at  which 
no  perforation  occurs  and  the  three  lowest  striking  velocities  at  which 
perforation  occurs.  These  data  are  then  averaged  to  give  a mean  value 
of  the  impact  velocity  which  is  the  mean  limit  velocity,  These 

data  may  be  or  may  not  be  mixed,  and  the  accuracy  of  the  estimate  will  be 
determined  by  its  range.  If  the  range  is  small  the  mean  limit  velocity 
can  bo  considered  reasonable  whether  the  response  is  mixed  or  not.  As 
the  range  increases  the  estimate  becomes  unreliable. 

B,  Ballistic  Limit  Velocity,  Calculation 


In  actual  experimental  tests  the  responses  are  generally  mixed  and 
only  a rough  estimate  of  the  limit  velocity  can  be  determined  by  the 
previous  method.  To  facilitate  the  analysis  of  sensitivity  data  Golub 
and  Grubbs^’®  of  this  laboratory  developed  a computational  procedure 
to  determine  the  ballistic  limit  velocity,  Vg^,  It  embodies  the  theory 

of  estimation  by  the  method  of  maximum  likelihood  developed  in  several 
uatistical  texts®.  The  procedure  accepts  sensitivity  data  in  any  form 
i'.nd  computes  the  maximum  likelihood  estimate  of  the  ballistic  limit 
velocity,  Vg^,  the  standard  deviation,  the  variance  and  covariance  of 

the  estimate,  and  the  confidence  intervals  for  selected  probability 
levels. 


VI.  ANALYSIS  OF  DATA 

Thirty  groups  of  data  were  selected  at  random  for  calculating 
critical  impact  velocity.  All  the  groups  had  Vj^-Vg  data,  and  their 

experimentally  determined  limit  velocities,  were  included.  Groups 

containing  only  quantal  response  were  not  considered  because  no  com- 
parisons could  be  made  of  their  ballistic  limits  by  another  method. 

®«/.  A,  Fevolii  "The  AaowKiaij  and  Heproduaibility  of  Seveml  Methods 
for  Obtaining  Ballistia  Limits  of  Amor"  Development  and  Proof  Servioesy 
Aberdeen  Proving  Qroundy  TB4-00bD,  AD~1Z00,  May  1S6?, 

^A,  OoUib  and  F,  E,  Orubbe,  "On  Estimating  Ballistia  Limit  and  Its 
Preoi8ion"y  Ballistias  Hesearan  Laboratories  Teahniaal  Note  Idi, 

Au  bll20y  March  1B60. 

S,  iiagan  and  y.  VisnaUy  "Analysis  of  Sensitivity  Data  Following 
A Nomal  DistriOution" y Analytical  Section  Report  70~A6-K3,  Materiel 
Test  DirectoratSy  October  1B70. 

®iV.  d,  Dixon  and  A,  M,  Moody  "A  Method  for  Obtaining  and  Analyzing 
Sensitivity  Data"y  Jourml  of  the  American  Statistical  Association, 

Vol.  46,  1D43. 
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The  mean  limit  velocity  Vg^g  was  calculated  for  all  the  groups  that  had 

no  mixed  response.  The  Vg^g  of  those  groups  that  had  loss  than  six  data 

points  is  indicated  by  a parenthesis  and  is  intended  only  as  an  approxi- 
mate value.  The  Vgy  was  calculated  for  all  the  groups  that  had  a 

mixed  response,  and  the  limit  velocity,  Vj^j^,wwas  computed  for  all  the 

data.  The  results  are  tabulated  in  Table  1. 


The  simplest  approach  yielded  a mean  limit  velocity,  Vg^g,  which 
approximated  that  of  the  limit  velocity,  Vj^.  Differences  in  the 


limit  velocities  ranged  from  1 to  28  meters/second  with  an  average 
difference  of  7 meters/second.  Since  all  the  values  of  Vg^g  were 


included  a mean  limit  velocity,  irrespective  of  the  number  of  data 
points  used,  will  compare  favorably  with  a limit  velocity.  The  only 
restriction  is  that  the  data  points  do  not  show  mixed  results.  If  the 
response  is  mixed  or  the  range  of  the  data  is  too  great  a mean  limit 
velocity  may  be  calculated  that  is  not  reasonable. 


The  determination  of  the  limit  velocity  by  the  Lambert  method 

yields  values  that  range  from  -36  to  +25  meters/second  with  an  average 
difference  of  7 raeters/second  with  80%  of  the  values  being  within 
4 meters/second  of  the  limit  velocity  Vj^. 

The  apparent  failure  of  the  Golub-Grubbs  procedure  to  compute  a 
VgQ  for  all  tests  is  attributed  to  several  restrictions  imposed  by  the 

procedure.  The  first  of  these  is  the  existence  of  a zone  of  mixed 
results  upon  which  the  estimate  of  Vg^  is  dependent.  A zone  of  mixed 

results  exists  when  the  highest  velocity,  called  BIGA,  at  which  the 
projectile  is  defeated  by  the  target  is  greater  than  the  lowest  velocity, 
called  SMALL,  at  which  the  target  is  defeated.  When  no  zone  exists, 
i.e.,  SMALL  > BIGA,  the  first  estimate  of  the  mean  is  zero,  and  other 


procedures  such  as  Vg^g,  or  must  be  employed. 


Ar.  alternate  use 


of  the  procedure  when  nc  zone  exists  is  to  provide  as  i.iput  an  initial 
value  for  the  mean  a. id  sigma,  and  to  vary  both  until  convergence  in  the 
iteration  process  is  obtained.  However  this  procedure  is  time  consum- 
ing and  it  is  easier  to  use  one  of  the  other  ‘.echniques.  When  BIGA 
equals  SMALL  its  velocity  is  the  true  Vg^.  Test  No.  4 deijonstrates 

this  condition.  When  a zone  of  mixed  results  exists  the  iteration 
process  which  calculates  the  best  estimate  is  dependent  on  .sample 
size.  If  BIGA  is  much  greater  than  SMALL  the  computation  will  tail. 

An  example  of  such  a condition  would  be  where  BIGA  is  the  largest  data 
point.  However,  a reduction  in  the  zone  of  mixed  results  through  an 
increase  in  the  number  of  samples  at  higher  striking  velocities  would 
insure  successful  computation  of  a Vg^. 
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till  nicoiui  rc^iritiion  i>  the  iaMnle  size.  The  c»timutc  of  V,,  is 

50 

itl.itivcly  iiidci'cnilfiil  of  sample  size,  but  the  estimate  of  standard 
‘h-vi.iiion  IS  not  only  biased  but  extremely  variable  for  small  sample 
•t:e.  iiul  nuj)-  provide  quest lonuble  conclusions  in  the  determination  of 
-t'MiiUence  intervals.  Kesults  for  a parametric  study,  using  'Jests  No. 
1 ' as  j typical  example,  demonstrates  the  fact  that  increasing  the 
•ampT  si..r  mil  provide  a solution  to  an  otherwise  meaningless 

. oiiipc  ! I ? 1 Ull. 
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Table  1.  Critical  Impact  Velocity  Results 


Test 

Sample 

V 

'^50S 

'^50 

No. 

Size 

m/s 

m/s 

m/s 

m/s 

1 

12 

832 

833 

807 

2 

17 

831 

843 

825 

3 

16 

860 

871 

858 

4 

12 

779 

780 

779 

5 

18 

847 

843 

849 

6 

12 

846 

892 

835 

7 

8 

890 

899 

888 

8 

7 

888 

896 

879 

9 

7 

1234 

1248 

1234 

10 

15 

732 

714 

720 

11 

17 

634 

(651)* 

625 

12 

20 

549 

(542) 

548 

13 

20 

569 

572 

575 

14 

20 

354 

358 

351 

15 

20 

335 

333 

334 

16 

17 

323 

322 

323 

17 

20 

284 

299 

293 

18 

7 

699 

(709) 

719 

19 

9 

858 

(867) 

858 

20 

8 

1113 

(1118) 

1149 

21 

4 

1283 

(1293) 

1299 

22 

10 

1452 

(1455) 

1450 

23 

15 

788 

810 

79? 

24 

8 

933 

(940) 

937 

25 

6 

1205 

(1201) 

1206 

26 

13 

1417 

1413 

1412 

27 

15 

1308 

1318 

1293 

28 

5 

990 

(993) 

994 

29 

5 

1140 

(1168) 

1158 

30 

10 

1296 

(1301) 

1296 

*(  ) means  Imit  velocity  determined  with  less  than  six  data  points. 
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VII.  CONCLUSION 


The  determination  of  the  method  to  obtain  an  accurate  critical 
impact  velocity  depends  upon  a careful  examination  of  the  data.  If 
no  zone  of  mixed  results  exists,  a V^^g,  Vj^  and  Vj^j^  can  be  computfd 


irrespective  of  sample  size.  If  a zone  of  mixed  results  exists,  the 
Golub-Grubbs  method  will  provide  a reasonable  value  of  Vg^  for  the 


critical  impact  velocity  for  a large  sample  size;  but  the  method  may 
fail  for  a small  sample  size  or  for  a very  large  zone  of  mixed  results. 
Furthermore,  in  conducting  any  experiments  where  the  results  desired  are 
in  the  form  of  sensitivity  data  it  is  advisable  to  use  some  statistical 
procedure  to  design  the  experiments  so  that  the  resulting  data  will  lie 
within  a restricted  range. 
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